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We observed the effect of radiation pressure on the angular sensing and 
control system of the Laser Interferometer Gravitational- Wave Observatory 



O 

O (LIGO) interferometer's core optics at LIGO Hanford Observatory. This is 

0> the first measurement of this effect in a complete gravitational wave interfer- 

^ ometer. Only one of the two angular modes survives with feedback control, 



since the other mode is suppressed when the control gain is sufficiently large. 
We developed a mathematical model to understand the physics of the system. 
This model matches well the dynamics that we observe, 
(c) 2009 Optical Society of America 

OCIS codes: 000.2780, 120.3180, 120.2230, 220.1140. 



1. Introduction 

Gravitational wave detection technology has made great strides in recent years, and more 
progress is expected in the near future. In the U.S., the Laser Interferometer Gravitational- 
Wave Observatory (LIGO) project has constructed several large interferometers, two at 
Hanford WA and another one at Livingston LA. [1] These initial LIGO interferometers were 
used to make observations of unprecedented sensitivity during the period November 2005 
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through September 2007, LIGO's fifth science run, during which one year's worth of triple 
coincident data were collected at an rms strain sensitivity of about 10~^^ [2-4]. The LIGO 
Scientific Collaboration (LSC) is now completing the analysis of these data. Operation of 
LIGO and analysis of gravitational wave data are closely coordinated with two European 
projects, GEO (a British-German collaboration with an 0.6 km interferometer near Han- 
nover [5]) and Virgo (a French-Italian collaboration with a 3 km interferometer near Pisa [6]). 

Striking as this technological progress is, it is widely expected that sensitivity improvement 
of about a factor of ten in strain amplitude will be required before observations of gravita- 
tional waves become routine. Fortunately, the technology to make the required improvement 
is in hand. For LIGO, those improvements will be embodied in Advanced LIGO, now under 
construction and expected to begin operations around 2015 [7]. In the meantime, one addi- 
tional data run will be carried out with an improved version of the initial LIGO hardware. 
This incremental improvement, called Enhanced LIGO, will run during 2009 and 2010. 

2. Gravitational wave interferometers 

Gravitational waves interact with matter by causing a distinctive strain (fractional separa- 
tion change) pattern in sets of test masses that are free to move. The distinctive pattern is 
that the sign of the strain is opposite in two perpendicular directions in the plane normal 
to the wave's propagation direction. For example, a wave propagating in the z direction 
causes equal and opposite strains in directions x and y. This means that a well designed 
Michelson interferometer can be a good sensor of gravitational waves. A properly aligned 
wave traveling in the vertical direction would cause simultaneous-but-opposite changes in 
the lengths of the two arms, thus yielding a phase difference between the light from the two 
arms that can be observed when the light recombines at the beamsplitter. 

Such a Michelson interferometer needs to include the following features: mirrors that are 
free to move in response to the wave (playing the role of the test masses) , isolation of those 
mirrors from other causes of motion other than gravitational waves (such as the vibration of 
the laboratory fioor) , large length (so that strain is converted to a sufficiently large relative 
motion between the mirrors) and sufficient phase sensitivity to allow the motions to be 
observed. The fact that strain sensitivity of 10~^^ has still not yielded a detection shows 
that these demands are severe. 

LIGO has addressed these challenges as follows. [1] (Other interferometers have solutions 
that are similar in most features but differ in details.) LIGO's interferometers have arms 
that are 4 kilometers in length. (The second of the two interferometers at LIGO Hanford 
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Observatory is only 2 kilometers long.) The 10 kg mirrors are suspended as carefully engi- 
neered pendulums, which are in turn connected to multi-stage seismic isolation systems. A 
high power Nd:YAG laser (with a power of 10 W) is used to provide sufficient optical power 
to reduce the phase noise due to photon shot noise. The simple Michelson interferometer 
configuration is augmented by several features that substantially enhance efficiency. Each 
arm is made up of a Fabry-Perot cavity held on resonance, thus amplifying the phase 
shift by a factor of JF/tt (the cavity's finesse T is about 200). The effective laser power 
is enhanced by an additional factor of about 50 by placing a partly transmitting mirror 
between the laser and the rest of the interferometer; the combination of this recycling mirror 
and the rest of the interferometer forms an additional Fabry-Perot cavity that is held on 
resonance. The whole system is enclosed in a vacuum system held at 10~* torr. (See the 
schematic diagram in Figure 7.) 

3. Role of radiation pressure 

Consideration of the LIGO design reveals that radiation pressure cannot be neglected in 
its operation. Light arrives at the recycling mirror with a power of 5 watts, but inboard 
of the recycling mirror the cfTcctivc power level hitting the interferometer's beamsplitter 
is 250 watts. Inside the resonant cavities that make up each of the long Michelson arms, 
the power level is 12 kilowatts. Thus, the radiation pressure on each of the arm mirrors 
is 10^^ ncwtons. Incident on a 10 kg mirror held in a pendulum suspension with resonant 
frequency of 0.5 Hz, this force would cause a displacement of 10^^ m. Of course, a Michelson 
interferometer with suspended mirrors only functions properly with a control system to hold 
the mirrors at the proper separations and alignment; this control system acts to counteract 
the radiation pressure force. A couphng to longitudinal motion of the mirrors has been 
experimentally demonstrated by Nergis Mavalvala et. al [8]. 

In 2006, Sidles and Sigg [9] pointed out that torques due to radiation pressure can cause 
an interesting problem in interferometers with suspended mirrors. They showed that the 
angular coupling between the mirrors due to radiation pressure is best understood by 
considering the dynamics of a cavity as a set of normal modes. At large enough radiation 
pressures, one of these modes can become unstable. Although LIGO detectors control not 
only longitudinal motion but also angular motion of the mirros, the servo bandwidth cannot 
be arbitrarily turned up due to noise coupling into the gravitational wave output. 

In this paper, we first review the theory developed by Sidles and Sigg, giving a heuristic 
understanding of the instabihty that they predicted. We point out that initial LIGO 
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operated successfully at power levels that would be expected to cause angular instability, if 
one made a naive application of their theory. In the next section, we present a measurement 
of the angular response to torque applied to a LIGO mirror, made under normal operating 
conditions (and at other light power levels as well.) The mirror response shows features 
different from those expected by the theory as presented in reference [9] . 

The difference from the simple theory is not surprising, since a control system acts on the 
mirrors of the interferometer. In order to derive the expected response, we modeled the 
combined mechanical-optical-control system using Matlab and Simulink [10]. Section 6 of 
the paper describes this model, and shows how our measurements agree with this more 
complete theory of the dynamics of the system. In particular, we can now understand the 
stability of initial LIGO at power levels at which instability was expected in an uncontrolled 
interferometer. Our measurements show that, inside the control system, the mirror's angular 
transfer function is indeed that of an unstable harmonic oscillator. 

In the final section of the paper, we investigate how much the light power can be increased 
before the present control system would be overwhelmed by the radiation pressure instability. 

4. Review of Sigg-Sidles instability 

Consider an optical resonator that consists of two suspended mirrors and laser light which 
resonates between the mirrors (Figure 1). The sketch shows two mirrors, which have very 
large radius of curvature (larger than the distance between two mirrors in LIGO detectors) ; 
the optical axis is the line connecting the centers of curvature of the mirrors. The two 
mirrors feel a natural restoring force that comes from the wire loop suspending them from 
the suspension system. 

When the cavity is resonant, radiation pressure from the laser in the cavity can become 
important. With no light, each mirror exhibits independent torsional oscillations. When the 
cavity is filled with fight, there are two coupled modes of the opto-mechanical system. (This 
description can be applied either to the yaw mode or the pitch mode of the mirrors in the 
cavity. The measurements described below were made in the yaw mode.) 

Suppose that the mirrors oscillate in such a way that the the sign of mirror angle ^1,^2 
is either (-I-, — ) or (— ,-|-). In this situation, radiation pressure works so as to enhance the 
original restoring force by pushing the mirrors back to the original position as shown in 
Figure 2, which means the mirrors experience a stronger restoring torque than that due to 
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Fig. 1. Two mirrors are coupled by radiation pressure caused by the power of 
the laser beam P, pushing at distance Xi and X2 away from the center line when 
the mirrors are tilted at angles of 62 respectively. Ii, I2 arc the moments of 
inertia of the mirrors. -Ri, i?2, L are the radii of curvature of the two mirrors 
and the distance between the mirrors. T is an external torque applied to one 
of the mirrors to study the system. The sketch can be viewed as either yaw 
motion or pitch motion. 



the mechanical restoring torque alone. Thus, the eigenfrequency of this mode will be higher 
than the original uncoupled pendulum natural frequency. 

On the other hand, suppose that both mirrors tilt in such a way that the sign of mirror 
angle ^1, 62 is either (+, +) or (— , — ). In that situation, the radiation pressure works against 
the original restoring force. In this case, if the tilt angle becomes larger, the beam spot 
will move as shown in Figure 3, which means that the net restoring force will be smaller 
than that due to the mechanical restoring force alone. When the power inside the optical 
resonator exceeds a critical value, the net restoring torque will become negative. 

In order to investigate the dynamics of this system, wc can apply an external torque to one of 
the mirrors and determine the transfer function of the mirror's angle to the applied torque. 
The normalized transfer function will be 
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stiff UJ+ > iuo 



Fig. 2. Stiff mode of angular motion in optical resonator. Radiation pressure 
works to enhance the mechanical restoring force. 
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where Uq is the mechanical resonant frequency of the torsion pendulum, P is the laser 
power inside the cavity, L is the length of the cavity, / is the moment of inertia of 
the mirror, and c is the speed of light. The factors gi and g2 are the g parameters of the 
cavity, defined by 511,2 = 1 — I//^i?i,2 where Ri and R2 are the radii of curvature of the mirrors. 



Here is how one can have an intuitive justification of the transfer function, which describes 
the process of shaking one of the mirrors and observing the response of that mirror. 
Since there are two angular modes of the coupled system, one will see resonances at the 
two different frequencies which correspond to the two pairs of poles uj±. Given the phase 
relations between these two modes, there will be a frequency at which the mirror docs not 
move at all. This frequency corresponds to a pair of zeros a;^. If we observe the other mirror. 
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soft UJ- < UOq 



Fig. 3. Soft mode of angular motion in optical resonator. Radiation pressure 
works against the restoring force from the wire, and if the power exceeds a 
critical value, then this mode is unstable. 

however, the pair of zeros does not show up. 

As P increases, uoj^ increases, while c<j_ and decrease. (Sec Figure 4.) Eventually, both uo"^ 
and uol become negative, and the corresponding poles and zeros become real. Pi and P2 in 
Figure 4 are the powers that give a;^ = and uj^ = 0, respectively. For initial LIGO, they 
are roughly Pi = 7.5 kW and P2 = 26.4 kW. A set of Bode diagrams are shown in Figure 5, 
and a pole map of (1) is shown in Figure 6. 

When P — 0, the two mirrors are uncoupled. The transfer function must be that for a simple 
torsion pendulum. When P is non-zero but sufficiently small, the two mirrors will be coupled 
and thus the Bode diagram will have two peaks associated with the two pairs of poles and 
one dip from the pair of zeros. In the s-plane, the two pairs of poles at s = ±iu!Q denoted 
by triangles in Figure 6 split into non-degenerate pairs, one moving toward the origin and 
one moving away. When P reaches Pi, the pair associated with a;_ reaches the origin, and 
for higher powers those poles move along the real axis. When P reaches P2, above which a 
pair of zeros become imaginary, the dip corresponding to cuz will disappear from the Bode 
diagram. The circles in Figure 6 are the poles at the operating laser power in initial LIGO. 
The existence of the pole in the right half plane indicates instability. At P > Pi, initial LIGO 
would be unstable in the absence of an angular control system. 
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Fig. 4. Angular frequencies of poles and zeros as a function cavity laser power, 
according to Equations (4.1) - (4.3). As P increases, cj+ increases, while uo^ 
and ujz decrease as P increases and become negative above Pi, P2 respectively. 



5. Measurement 

To study the radiation pressure effect, we made a set of measurements on one of the arm 
cavities in the 4 km interferometer at LIGO Hanford Observatory. We wanted to study the 
regime of P > Pi, as described above, so all control loops (including angular controls) were 
engaged during the measurement. The angular sensing and control system [11-14] has ten 
degrees of freedom. They are the pitch and yaw motion of the five core optics (the two pairs 
of mirrors: (ITMX, ETMX) and (ITMY, ETMY) in the long Fabry-Perot cavities, and the 
recycling mirror (RM) located upstream of the Michelson interferometer, see Figure 7.) 
These degrees of freedom are measured by the quadrant photodetectors called wavefront 
sensors, [15, 16] and are controlled by electromagnetic actuators attached to the backs of 
the mirrors. We reduce the ten degees of freedom to two by looking at only two wavefront 
sensors which are sensitive to the differential degrees of freedom, (i.e., A.9etm — O2 — O4, 
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(C) Pi < P < P2 



(d) P2 < P 



Fig. 5. Bode plot of the uncontrolled two-mirror system at various powers, (a) 
P = 0, (b) < P < Pi, (c) Pi < P < Pa, (d) Pa < P. Note that the bottom 
two plots correspond to the unstable cases where external servo is necessary 
to keep the system stable. 



and AOjTM = - ^3) 



Besides these wavefront sensors for global angular sensing, there are sensors called optical 
levers [17], which work only locally. An optical lever consists of a diode laser and a quadrant 
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Fig. 6. Poles in the s-plane at various powers, for the uncontrolled system. 
Triangles on the imaginary axis at P = split into two directions: uj+ moves 
away from the origin while c<j_ moves toward the origin. After the pair for a;_ 
hits the origin at P = Pi, the poles move apart along the real axis. 



photodiode which senses the position of the laser spot reflected by the mirror, and thus 
allows a measurement of the angle of the mirror. We use the optical lever signal to monitor 
the response of the mirrors to angular excitation; this enables us to measure the transfer 
function H{s) that we investigated in the previous section. 

In order to start the measurement, we flrst ensure that the interferometer is in its ordinary 
operating condition, where the two arm cavities and the recycling cavity are resonant, and 
the complete control system is engaged. Then, we inject an excitation signal into one of 
wavefront sensors control path that results in dithering the end test masses differentially. 
During the excitation, we monitor the optical lever signal of the end mirror in the x-arm 
and a signal which goes into the driver of the actuator attached to the mirror. Since the 
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Fig. 7. Block diagram of our measurement scheme. An excitation signal is 
injected into the WFSl path that dithers ETMX and ETMY differentially, 
i.e., A9etm = 02 — 64. Both WFSl and WFS2B are sensitive to the differential 
DOF {A9etm, AOitm), but we focus on the contribution to the x-arm. The 
symbols next to each core optic show the rotation axes for pitch and yaw tilts. 
A positive tilt angle follows a right-handed rotation about the axis. 



driver signal is proportional to the torque produced by the actuator, we are able to monitor 
the transfer function of the mirror angle D to torque C applied to the same mirror up to 
some overall gain. The transfer function is the direct analog of the one we introduced in 
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the previous section, i.e., the response of the ETMX mirror angle to the torque apphed 
to the mirror. In addition, we record signals from two points located just upstream A and 
downstream B of the excitation point in the wavefront sensors control loop. This enables us 
to calculate the open loop transfer function of the control loop. We here call the two transfer 
functions oplev{s) and olg{s) respectively. 

oplev{s) = ^ (4) 

olg{s) = ^ (5) 

We made the measurement at three different input laser powers (0.8 W, 4.0 W, and 6.8 
W.) [18]. In Figure 8, the open symbols show the results of our measurements of the oplev{s) 
transfer function, while Figure 9 shows our measurements of olg{s). (The solid lines are the 
predictions of our mathematical model including the effect of control loops, described in the 
next section.) 

There are several surprising features in Figure 8. Firstly, note that there is only one peak 
in the transfer function at each laser power, whose frequency decreases as the laser power 
goes up. This is something that we did not expect based on the naive model (i.e., the model 
without any control loops). Note also that, at the highest power level, a 180 degree phase lead 
was observed. This was also surprising since it would usually indicate the system's instability. 



6. Mathematical model 

These observations led us to construct a more complete model of the system that explicitly 
includes the dynamics of the angular control system. Figure 10 shows the block diagram. 
The most important parts of the control loop are the sensing matrix (labeled as a, 6, c, d) 
and the control matrix {B,C,D). Each wavefront sensor measures a combination of the 
angle of the two mirrors (input and end mirrors); the signals are fed back to each mirror 
to minimize the deflection angle. The sensing matrix is measured separately in advance 
by shaking each mirror; the control matrix is determined by inverting the sensing matrix. 
Besides the wavefront sensors for global angular sensing, there is also a local loop called 
optical lever for each mirror, and a compensation loop to make the local loop invisible to 
the wavefront sensors. All filter banks associated with the two wavefront sensors arc built 
into the model. The solid curves in Figures 8 and Figure 9 show the transfer functions 
calculated by our models. The different colors in the plots denote different cavity laser 
power levels in the model. They are blue for 1.6 kW, green for P = 7.0 kW, and red 
for P = 12.5 kW. The agreement between the measurement and the model is very good. 
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Frequency [Hz] 



Fig. 8. Bode plot of the optical lever transfer function, oplev{s). The input 
laser power levels are indicated as follows: 0.8 W -.circles, 4.0 W -.triangles, and 
6.8 W -.squares. The solid lines are the predictions of our mathematical model 
including the effect of control loops. 

The only adjustable parameter in the model is the cavity laser power, (related to the 
input power by an imperfectly determined factor.) Knowing that the full interferometer has 
many more degrees of freedom than are modeled here, such agreement is especially gratifying. 

There are several features worth noting in these transfer functions. Firstly, we can now 
understand why there is only one peak observed in the optical lever transfer function. This 
resonance is associated with a pair of poles s = ±zu;_, but the control loops make the 
transfer function much more complicated than (4.1). The reason why we only have one 
visible resonance is that with high loop gain as in initial LIGO, the pair of zeros s = ±iujz 



13 




Fig. 9. Bode plot of the open loop transfer function, olg{s). 



and the other pair of poles s = ±2u;+ asymptotically approach the same value and cancel 
each other. 

Secondly, we can now understand why the phase of the optical lever transfer function at 
6.8 W has a 180 degree phase lead. As noted above, under the influence of a strong control 
loop, the pair of poles associated with the stable mode cancels the pair of zeros and only the 
pair of unstable poles survives. In other words, when the control loop strongly suppresses 
the motion of the test masses, the in-loop transfer function oplev{s) looks as if it were a 
single mode unstable oscillator. This is clearly observed in both our measurement and the 
model as the phase advancement in the plot. This is evidence that the radiation pressure 
effect on optical resonators, predicted by Sigg and Sidles, really exists in the interferometer 
and that the system successfully operates in the regime where the instability would occur 
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Mechanical Part 




Fig. 10. Block diagram of interferometer arm with angular feedback con- 
trol, a, b, c, d and B, D, E are the sensing and control matrix, respectively. 
Rectangular blocks stand for filter banks. The SUS filters are ones related 
to the suspension system. The factor Kp stands for the optical parameter 
2PL/ [cl{l-g,g2)]. 

in the absence of control. 

Thirdly, the frequency associated with the peak shifts toward lower frequencies as laser 
power goes up. This is understood as a radiation pressure effect, since the peak we observed 
corresponds to a pair of poles s = ±iuj_ and uj_ decreases as shown in Figure 4. The 
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frequency does not go down as rapidly with increasing power as in the naive model, since 
the control system enhances the effective restoring force on the ITM. 

Finally, we note that the radiation pressure effect is also observed in the phase of the open 
loop transfer function olg{s). 

7. Stability 

It is very important to know how much more power can be handled with the current control 
scheme for the upcoming higher power operation of the LIGO interferometers. In particular, 
an upgrade to initial LIGO (called Enhanced LIGO) is now being developed. It will have 
laser power roughly three times greater than initial LIGO, and will have an almost identical 
control system to that in initial LIGO. 

When we discussed the stability of the system without control, a torque was applied to the 
system from 'outside'. (We called it T.) However, with the control system in place, T is no 
longer outside, but is included within the system. Therefore, in order to judge the stability 
of the whole system, we need to consider another torque T', which is located truly outside 
of the whole system. (See Figure 10.) We will investigate the response of the system to such 
an outside torque 

H*{s) = ^. (6) 

Figure 11 shows how the system poles move in the s- plane as a function of the cavity laser 
power. We immediately notice the difference between Figure 11 and Figure 6. The circles, 
which indicate the systems poles at the full laser power level for initial LIGO, are now all 
in the left half plane, as expected. According to our model, the system will remain stable 
until the cavity laser power reaches roughly P/Pligo = 8.5. 

This analysis was based on the degrees of freedom associated with differential motion 
between the mirrors in the interferometer's two arms. The common mode degrees of freedom 
have lower unity gain frequencies in LIGO. This means that the full initial LIGO system 
would probably not remain stable at powers as high as are predicted in our analysis. A more 
complete analysis including the common mode degrees of freedom will be left for future 
work. We also note that a new approach on this issue is being tested during the Enhanced 
LIGO commissioning process [19]. It uses a control scheme built around the stiff and soft 
modal basis. 
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Other workers in the field have also made studies of the radiation-pressure-induced angu- 
lar instability in optical cavities. Sakata has demonstrated the effect in a specially-built 
apparatus at the National Astronomical Observatory of Japan [20]. 
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Fig. 11. Poles of the system with feedback control, with laser power from zero 
to 110 kW. The power increment between points is 2.5 kW. Triangles and 
circles are for P = 0, P = Pligo — 12.5 kW, respectively. 



8. Summary 

We observed the effect of radiation pressure on the angular control of the LIGO core optics 
at the 4 km interferometer at LIGO Hanford Observatory. This is the first measurement of 
this effect performed on a full gravitational wave interferometer. Only one of two angular 
modes survives with feedback control since the other mode is suppressed when the control 
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gain is large enough. A mathematical model was developed to understand the physics. It 
indicates that the system will remain stable at substantially higher power levels. 
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